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Abstract —Wind is becoming an increasingly significant 
source of energy in modern power generation. Amongst ex¬ 
isting technologies, Variable Speed Wind Turbines (VSWT) 
equipped with Double Fed Induction Generators (DFIG) 
is widely deployed. Consequently, power systems are now 
experiencing newer power flow patterns and operating con¬ 
ditions. This paper investigates the impact of a DFIG based 
Wind Farm (WF) on the stability of electromechanical 
oscillations. This is achieved by performing modal analysis 
to evaluate the stability of a two-area power network 
when subjected to different wind penetration levels and 
different geographical installed locations. The approach via 
eigenvalues analysis involves the design of voltage and Sup¬ 
plementary Damping Controllers (SDCs) that contribute to 
network damping. The effect of Power System Stabilizer 
(PSS) is also examined for several network conditions. 
Simulations demonstrate a damping improvement up to 
933% when the control systems are activated and the system 
operates with 25% wind integration. 

I. Introduction 

Today, DFIGs are the most popular wind power genera¬ 
tors in the wind power industry. The use of this machine 
has certain advantages, such as variable speed, active- 
reactive power control (by using converters) and the 
economic cost Q. The layout of a DFIG Wind Turbine 
(WT) is shown in Fig. [I] 0. Between the rotor and the 
grid there are two linked converters in a back to back 
pattern. The only requirement for the power electronic 
converters is the percentage that must be rated in order 
to handle a fraction of the total power (rotor power). In 
most cases, this percentage is about 30% of the nominal 
generator power 0 . The two converters are connected 
through a dc-link capacitor which, as an energy storage 
element, maintains small voltage variations in the dc-link. 
The control of the torque (or the speed) of the DFIG as 
well as the power factor at the stator terminals is achieved 
with the Rotor Side Converter (RSC). The main objective 
for the Grid Side Converter (GSC) is keeping the dc-link 
voltage constant despite of the magnitude and direction of 
the rotor power. The GSC operates at the grid frequency, 


primarily to produce (leading) or to consume (lagging) 
reactive power. Typically, a transformer is connected in 
the middle of the GSC, or the stator, and the grid. The 
working frequency of RSC depends on the wind speed 
0 . 

Fig. [2] shows the equivalent model of the turbine 
according to Parks equations, referenced to the stator. The 
circuit contains the speed voltage vj'ijj and the stator and 
rotor components m. 

The authors in 0 discuss how wind power and the 
type of the WT generator change the damping and the 
frequency of eigenvalues. In particular, the study shows 
that wind power tends to improve the damping of inter¬ 
area oscillations while the impact on intra-area is not 
significant. In 01 it is shown that New Zealand systems 
damping performance is not significantly affected by 
increasing the wind penetration level. An efficient way 
to increase the damping of electromechanical oscillations 
using an auxiliary PSS loop installed in DFIG WT 
generators is presented in 0. Fernandez et al. 0 add 
controllers in active and reactive power control loops 
of DFIG to enhance the damping of oscillations. In 
addition, synchronous generators’ replacement with DFIG 
based WF in Q, it is shown that does not degrade the 
two area system modal characteristics. Nevertheless, WF 
improves the stability of the inter-area mode of the power 
system. A damping control method for a STATCOM in 
a series compensated WF for mitigating subsynchronous 
resonance and damping oscillations is analyzed in 0. In 
|9l non-linear controllers are implemented for improving 
the damping of intra-area and inter-area oscillations. 

The above references verify the impact of wind genera¬ 
tions on electromechanical oscillations. In most cases the 
goal is to design controllers that can improve network 
damping. However, the control methods presented in 
literature can be implemented only for certain systems 
and do not have universal usage. Moreover, the WFs 
consist hundreds of WT generators distributed on vast 
areas. Consequently, it is perplexing to select the feedback 
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Fig. 1. DFIG based WT generator (2). 



Fig. 2. DFIG equivalent circuit (T]|. 


signals. Also, conventional PSS methods are practically 
difficult to implement into the WT generators since their 
capacity is small and their number is greater compared 
to the synchronous machines. These differences between 
conventional power plants and WFs necessitate the design 
of a controller for DFIG based WFs which has simple 
theory and can be easily implemented. 

The layout of this paper is as follows: It begins by 
explaining the concepts of RSC, GSC and the dc-link of 
back-to-back converter. Section [TIT] describes the theory 
behind modal analysis. In Section [IV] the voltage and 
damping control strategy is presented. Section |Y] presents 
the test system and the study scenarios. Finally in Section 


VI simulation results are presented and comments on the 


results of the system performance are addressed. 


II. Description of the back-to-back converter 

This section discusses the arrangement of the back- 
to-back converter. It presents the rotor-side and grid-side 
current control reference frames and explains the utility 
of the dc-link. 


A. Rotor side converter 

The RSC applies the voltage to the rotor windings of 
the DFIG and its purpose is to control the rotor currents. 


Specifically, it controls the developed torque at the shaft 
of the machine as far as the rotor flux position is oriented 
in an optimal way with respect to the stator flux. 

On the rotor perspective, the current is disintegrated in 
two mechanisms: in phase with the stator flux linkage 

and in quadrature. The electromagnetic rotation and 
the machine excitation are decoupled and measured by 
means of i$ and correspondingly. Fig. [ 3 ] (a) shows 
the association between the rotor current phasor I r , the 
stator flux linkage ^ s , the current components used by 
the control system 1 $ and i^} 9 a synchronous rotating 
reference frame rf, q which formulates an angle (ip) with 
the stator flux linkage, and the current components in such 
reference frame id r and i qr lITOl . 

The RSC regulates the WT output power and voltage 
(or reactive power) at the machine’s stator terminals 
by using a torque controller such that the power will 
follow a turbine power-speed characteristic (tracking the 
maximum power point). The actual power from the gen¬ 
erator terminals ( P ou t X added to the total power losses, 
is compared with the reference power (P re f) obtained 
from the WT characteristic. In most cases, a Proportional- 
Integral (PI) regulator is used at the outer control loop to 
reduce the error (power/rotor speed error). The output of 
this regulator is iqfi i.e. the reference rotor current that 
must be injected into rotor winding by the RSC. 

B. Grid side converter 

The GSC is designed in such way to regulate the 
capacitor voltage of the dc bus. In addition, it can generate 
or absorb reactive power to provide voltage support 
whenever necessary. 

On the grid side, the current is also disintegrated in 
two mechanisms: in phase with the stator voltage 

and iqa in quadrature. The active and reactive power 
exchanged with the grid are also decoupled and controlled 
by means of and i q £ J, respectively. In this kind of 
control systems, the active power output (. P ou t ) of the 
GSC is used to manage the dc link voltage level. 

Fig. [3] (b) shows the connection between the stator 
voltage u s , the GSC current I a , the control system 
current components i^) and i q £ J, a synchronous rotating 
reference frame d, q which formulates an angle (s) with 
the stator voltage, and the current elements in this kind 
of reference frame %d a and i qa 11108 . 

C. DC-Link 

The dc-link capacitor variations of the voltage are 
represented by the dc-link model which takes as input the 
power. Eq. [I] shows the energy stored in the dc capacitor 

m. 
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Fig. 3. (a) Rotor-side current control reference frame, (b) Grid-side 
current control reference frame fTol . 


w dc = J P dc dt = \-C-Vl (1) 

where Wd c is the stored energy and Pd c is the de¬ 
link input power. Eq. [2] presents the energy and voltage 
derivatives: 

dVdc _ Pdc 
dt C • Vdc 

dW dc _ p U 

dt dc 

The P dc is calculated as Pd c = Pin ~ Pc • Pin is the RSC 
input power and P c is the GSC output power. The voltage 
at the dc-link is constant when Pd c = 0 and varies with 
respect to Pdc • 

The GSC controller monitors the dc-link voltage which 
essentially determines the exported power by the GSC. 
For example, if the dc-link voltage increases then the 
GSC can export more real power by increasing the load 
angle such that the voltage will reduce back to its nominal 
value. In simple terms, the dc-link voltage represents 
the power flow balance between generated and exported 
energy in the rotor side. In the case where the dc-link 
capacitor input and output power do not match, then the 
dc-link voltage will change accordingly lUTTl . 

III. Modal Analysis 

Modal analysis of angle stability is concerned with 
the characteristic modes determination of a system mode 
linearized about its Stable or Unstable Equilibrium Points 
(SEP, UEP). It primarily comprises computing eigenval¬ 
ues, eigenvectors and participation factors. Particularly, 
analysis of these quantities uncovers the involvement 
of those machines that are losing synchronism, as for 
example in an inter-area mode oscillation lfl2) . Modal 
identification relates to the perseverance of characteristic 
modes from large-disturbance dynamic behavior, that is 
gained through transient stability simulations. 


The non linear model that defines the dynamics of a 
power system, includes a set of differential and algebraic 
equations (Eq. [3] and [4| lH3l : 

x = f(x,y) (3) 

0 = g(x,y) (4) 

where x and y are the vectors of state variables and 
algebraic variables, correspondingly. 

A conventional modal analysis is implemented by 
linearizing the system equations in Eq. [3] about an equi¬ 
librium point, x 0 , defined by: 

Xo = f{x a ) = o (5) 

Small disturbances around x Q can then be defined by: 

Ax = A • Ax (6) 

where A is the system state matrix defined by: 

A = ( A ij) = (^rM=x 0 ) (7) 

The characteristic modes of Eq. [6] are generally formed 

as: 

Vi ■ e Xit (8) 

where pi is the characteristic vector (or right eigenvector) 
and is the corresponding eigenvalue. Eigenvalues are 
solutions of the characteristic equation: 

det(A - XI) = 0 (9) 

Lyapunov stability theory states that if all the eigenval¬ 
ues of matrix A have negative real parts, then the equilib¬ 
rium point is small signal stable lH3l . If the eigenvalues 
are complex and located near the imaginary axis of the 
complex plane, then are characterized as critical. For an 
eigenvalue X = a ± juo, the frequency of the oscillation 
is given by uj/ 2n and the damping ratio ( is given by: 

( = —cr/V a 2 + uj 2 (10) 

The mode shape related to each critical eigenvalue can 
be used to identify the oscillation types and therefore the 
generators’ location in the power system |[T3l . In a power 
system with multiple generators, the dominant electrome¬ 
chanical oscillation mode often present two clusters rel¬ 
ative swing. Furthermore, when the electromechanical 
system is being low-frequency oscillation, the increment 
of angular speed and power angle of two generators have 
the characteristics of sinusoidal oscillation with same 
frequency and reversed phase. 
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Fig. 4. Voltage Controller. 


IV. Voltage and damping control systems 

A. Voltage Control 

The high penetration level of wind generations affects 
the power flows and hence the node voltages fl4l . Due to 
the fact that WFs replace conventional generators, WFs 
must contribute to voltage control. As a result, technical 
regulations for wind power generation are growing at the 
same rate at its installation El. The basic approach for 
the voltage control is to be centralized at the Point of 
Common Coupling (PCC) that is the voltage bus to which 
the WF is connected. Fig.[4]illustrates the configuration of 
the voltage controller. Voltage is transformed from three- 
phase to direct and quadrature (dq) components. The 
reference voltage that we want to maintain to the voltage 
bus is added to (dq) voltage and through a PI controller 
we get the direct axes current. This current controls the 
dq voltage of the RSC and hence we are able to regulate 
the voltage to the desired level. 

B. Supplementary damping control 

The unbalanced rotor torque of synchronous generator 
in power systems has the effect to cause electrome¬ 
chanical oscillations. An effective solution is adding a 
supplementary damping signal in the active power control 
loop. The configuration of the SDC system with the tradi¬ 
tional lead-lag compensation structure which contributes 
to improve the damping ratio of machines torque is shown 
in Fig. [5j It contains a washout filter, a gain block and 
two lead-lag blocks. The transfer function of the SDC 
system is presented in Eq. |TT] 


V pss (s) = [K] • 


sT w 


sTi n + 1 


sT 2 n + 1 

_sT w +1_ 


_ sTid + 1 _ 


_ sT 2 d + 1 _ 


(ID 


where K is the SDC system gain, T w is the washout time 
constant and T\ n , Tid, T 2n , T 2 d are the time constants of 
the lead-lag blocks. The washout time constant prevents 
steady changes in speed from modifying the field voltage, 
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Fig. 5. Supplementary Damping Controller. 


whereas the time constants of lead-lag blocks provide the 
necessary phase compensation. 

The design of the SDC system for WTs is based on the 
PSSs used for conventional generators. The input can be 
any signal affected by the oscillation. This fact implies 
the selection of the PCC as measurement point in order 
to avoid the filtering effect introduced by the transformer 
which is connected between the grid and the WF. In this 
study, the rotation rate deviation A uo of WF at the PCC is 
used as input signal. The amplification block affects the 
amplitude frequency characteristic. The role of washout 
block is filtering the direct current component. The lead- 
lag block is used to compensate the phase lag caused 
by the converter. The output signal can be any variable 
capable of varying the power delivered to the grid. For our 
case the excitation voltage is used as the output signal. 

The active power transfer function of the excitation 
controls with SDC system in the current regulator of the 


RSC control system is given by Eq. 12 


G{s) = ^-TT = C(sl - A)~ l B 
K ’ A u(s) v ' 


( 12 ) 


V. Test system and study scenarios 
A. Test system 

To evaluate the impact of a DFIG based WF on 
the electromechanical oscillations, the two-area power 
network shown in Fig. © El is examined. The sys¬ 
tem consists two identical areas connected through a 
relatively weak tie. Each area includes two synchronous 
generators with equal power outputs and loads which are 
modeled as constant impedance (full symmetry). Three 
electromechanical oscillations are present in this system. 
Two intra-area modes (i.e. generators’ oscillations in a 
certain area of the network against each other), one in 
each area, and one inter-area low frequency mode (i.e. 
generators’ oscillations in a certain area of the network 
against generators in another area of the network). 

The parameters of the study include different wind 
penetration levels and different geographical installed 
locations. The examined level of wind generation is up 
to 35% of the parallel synchronous machine capability. 
Locations constitute the connection of a WF next to the 
synchronous generator G4 and next to load LI as shown 
in Fig. [ 7 ] and Fig. [8} respectively. 
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Fig. 6. Kundur’s two-area system. 



Fig. 7. Kundur’s two-area system integrated with a WF next to the 
synchronous generator G4. 


Fig. 8. Kundur’s two-area system integrated with a WF next to the 
load LI. 


TABLE I 

Results of two-area four machine system 


Case 1 

Eigenvalue 

Damping 

Freq. 

(Hz) 

Oscillation 

mode 


0.108 + 4.027i 

-0.027 

0.641 

Inter-area 

No PSS 

-0.677 ± 7.048i 

0.096 

1.122 

Intra-areal 


-0.669 ± 7.269i 

0.092 

1.157 

Intra-area2 


-0.568 ± 4.055i 

0.139 

0.645 

Inter-area 

PSS 

-1.983 ± 8.276i 

0.233 

1.317 

Intra-areal 


-3.011 ± 9.764i 

0.295 

1.554 

Intra-area2 


B. Study scenarios 

In the case studies, four scenarios are considered. 
For each case study, the system is simulated with and 
without PSS for the synchronous conventional generators. 
Furthermore, a 3-phase fault takes place at the tie-line 
between the east and west area of the system at 1.0s and 
it vanishes 0.2s later: 

Case 1: Without WF connected to the system (Fig. [6]). 
Case 2: With WF integrated to the system next to the 
synchronous generator G4 (Fig. [ 7 ]). In order to distinctly 
highlight the fundamental characteristics of WF, no con¬ 
trols are applied to the WF. 

Case 3: With WF integrated to the system next to the 
load LI (Fig. [§. Like the previous case, no controls are 
applied to the WF. 

Case 4: With WF included to the system next to the load 
LI (Fig. [8]). In this case, voltage controller and SDC are 
integrated to the system. 


VI. Simulation Results 


The validation of the presented methods is performed 
by extracting the eigenvalues of the linearized state-space 
model using MATLAB/Simulink. The results regarding 
inter-area and intra-area oscillations are shown in Tables 
mm and [TV] for the Cases 1, 2, 3 and 4 respectively. 

From Tables id 0® and [TV] we observe that when 
the power system is integrated with WF, the damping 
of the system turns weaker without PSS control into 
the synchronous generators. However, when the PSS 
control is activated the inter-area oscillations are damped. 


Comparing this observation when voltage and SDCs are 
added, the system becomes more stable and the damping 
is improved. Particularly, in cases 1 to 3 i.e. without 
integrating voltage controllers and SDCs, the inter-area 
oscillations have a negative damping contribution when 
PSS is deactivated. Nevertheless, in case 4, even without 
PSS activated, the controllers provide positive damping 
to the system at the inter-tie frequency. The comparison 
between the test cases when PSS is used shows a 62% 
improvement when WF is attached next to LI and all the 
controllers are switched-on. Even in the situation where 
voltage and damping controllers are deactivated, there is 
a 18% improvement to the damping ratio with 35% wind 
integration next to G4. In order to show the great effect 
of the controllers, we choose the comparison between the 
base system (Fig. [6]) and the system presented in Fig. [8] 
since the latter has slightly worse inter-area oscillation 
than the system shown in Fig. [ 7 ] The improvement by 
using the controllers is around 49%. Finally, the overall 
improvement between the base study and the study where 
a) all the controllers are activated and b ) 25% wind 
generation is integrated to the system, is around 933% 
at 0.64Hz. 


VII. Conclusions 

In this paper, the impact of DFIG based WF on the 
small-signal stability performance is investigated. For 
the conditions considered, the WF does not degrade 
the two-area system modal characteristics when the PSS 
controller is activated into the conventional machines. 
On the contrary, PSS control improves the damping of 
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TABLE II 

Results of two-area four machine system with WF 
ATTACHED NEXT TO G4 


Case 2: 
No PSS 

Eigenvalue 

Damping 

Freq. 

(Hz) 

Oscillation 

mode 

10%WG 

0.094 ± 4.168i 
-0.679 ± 7.048i 
-0.625 ± 7.436i 

-0.023 

0.096 

0.084 

0.663 

1.122 

1.183 

Inter-area 

Intra-areal 

Intra-area2 

25%WG 

0.089 + 4.320i 
-0.680 ± 7.048i 
-0.483 + 7.720i 

-0.021 

0.096 

0.062 

0.688 

1.122 

1.229 

Inter-area 

Intra-areal 

Intra-area2 

35%WG 

0.084 ± 4.405i 
-0.328 + 7.934i 
-0.680 ± 7.049i 

-0.019 

0.041 

0.096 

0.701 

1.263 

1.122 

Inter-area 

Intra-areal 

Intra-area2 

Case 2: 
PSS 

Eigenvalue 

Damping 

Freq. 

(Hz) 

Oscillation 

mode 

10%WG 

-0.617 ± 4.221i 
-1.802 ± 8.394i 
-3.011 ± 9.763i 

0.145 

0.210 

0.295 

0.672 

1.336 

1.554 

Inter-area 

Intra-areal 

Intra-area2 

25%WG 

_0.694 ± 4.410i 
-1.463 ± 8.577i 
-3.011 ± 9.761i 

0.155 

0.168 

0.295 

0.702 

1.365 

1.554 

Inter-area 

Intra-areal 

Intra-area2 

35%WG 

-0.751 ± 4.522i 
-1.179 ± 8.703i 
-3.012 ± 9.760i 

0.164 

0.134 

0.295 

0.720 

1.385 

1.553 

Inter-area 

Intra-areal 

Intra-area2 


TABLE III 

Results of two-area four machine system with WF 

ATTACHED NEXT TO LI 


Case 3: 
No PSS 

Eigenvalue 

Damping 

Freq. 

(Hz) 

Oscillation 

mode 

10% WG 

0.141 ± 4.127i 
-0.604 + 7.180i 
-0.665 ± 7.278i 

-0.034 

0.084 

0.091 

0.657 

1.143 

1.158 

Inter-area 

Intra-areal 

Intra-area2 

25%WG 

0.172 + 4.156i 
-0.331 ± 7.104i 
-0.762 + 7.115i 

-0.041 

0.047 

0.106 

0.661 

1.131 

1.132 

Inter-area 

Intra-areal 

Intra-area2 

35%WG 

0.192 ± 4.173i 
-0.131 ± 7.058i 
-0.834 ± 7.01 li 

-0.046 

0.019 

0.118 

0.664 

1.123 

1.116 

Inter-area 

Intra-areal 

Intra-area2 

Case 3: 
PSS 

Eigenvalue 

Damping 

Freq. 

(Hz) 

Oscillation 

mode 

10% WG 

-0.603 ± 4.156i 
-1.984 + 8.276i 
-2.840 + 10.057i 

0.144 

0.233 

0.272 

0.662 

1.317 

1.601 

Inter-area 

Intra-areal 
Intra-area2 

25%WG 

-0.638 + 4.172i 
-2.123 ± 8.176i 
-2.633 + 9.961i 

0.151 

0.251 

0.256 

0.664 

1.301 

1.585 

Inter-area 

Intra-areal 
Intra-area2 

35%WG 

-0.660 + 4.18 li 
-2.221 ± 8.111i 
-2.473 ± 9.893i 

0.156 

0.264 

0.243 

0.665 

1.291 

1.575 

Inter-area 

Intra-areal 
Intra- area2 


inter-area modes when WF penetration level is increasing. 
The results obtained also show that when the voltage 
and SDCs are added to the DFIG the system inter-area 


damping ratio improves dramatically. 

TABLE IV 

Results of two-area four machine system with WF 
ATTACHED NEXT TO LI WITH CONTROLLERS AND 25% WIND 

Generation 


Case 4 

Eigenvalue 

Damping 

Freq. 

(Hz) 

Oscillation 

mode 


-0.060 ± 4.223i 

0.014 

0.672 

Inter-area 

No PSS 

-0.351 + 7.029i 

0.050 

1.119 

Intra-areal 


-0.760 ± 7.117i 

0.106 

1.133 

Intra-area2 


-0.933 + 4.040i 

0.225 

0.643 

Inter-area 

PSS 

-2.121 ± 8.173i 

0.251 

1.301 

Intra-areal 


-2.709 ± 9.828i 

0.266 

1.564 

Intra-area2 
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